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Abstract 

The changes in the electronic properties of the substitutional^ disordered MgC(Nii_ x T x ) 3 
(T=Fe, Co or Cm) alloys are studied using the atomic sphere formulation of the 
Korringa-Kohn-Rostoker coherent-potential approximation method (KKR-ASA CPA), 
while the effects of incipient magnetism in these alloys are studied phenomenolog- 
ically using Ginzburg-Landau coefficients in conjunction with fixed-spin moment 
method. We find that the disordered M gC (N ii_ x T x ) 3 alloys have a small magnetic 
moment localized at Fe and Co sites for low concentrations. The overestimation of 
the calculated magnetic moment is likely to be due to the limitations of the local- 
density approximation used in the present study. However, the calculated Ginzburg- 
Landau coefficients clearly show that the disordered MgC{Ni\^ x T x ) 3 alloys remain 
paramagnetic. At expanded volumes, we also find the possibility of a ferromagnetic 
state for MgC(Ni .Q 5 Feo m ) 3 and MgC(Nio, go Co . w ) 3 . 



Incipient magnetism akin to spin-fluctuations is anticipated [1,2] in the 8K 
superconductor MgCNi 3 [3]. A strong van Hove singularity just below the 
Fermi energy Ep [4,5] in MgCNi 3 suggests that the material may be close 
to a magnetic instability. Using a rigid band picture, it is estimated that 
approximately 0.5 holes would make MgCNi 3 magnetic [4,6]. However, in 
experiments the hole dopings accomplished via Fe or Co substitutions in the 
Ni sub-lattice of MgCNi 3 do not lead to magnetic solutions [7,8,9,10,11]. 

Model Hamiltonian- based calculations [12] suggest that dilute impurities of 
3<i transition-metals would develop a small magnetic moment in MgCNi 3 . 
Other calculations [13,8,14,15,16], based either on supercells or on ordered 
MgC(Nii_ x T x ) 3 alloys with T=Fe or Co, infer that in low concentrations the 
impurities would remain paramagnetic. The results of the Model Hamiltonian- 
based calculations are inconsistent with the experiments, at least for Fe and 
Co. The results of other calculations using a rigid-band-like picture, although 
consistent with the experiments in so far as the lack of magnetism is concerned, 
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do not provide a reliable and a complete description of substitutional^ disor- 
dered MgC(Nii_ x T x ) 3 alloys. The present study is intended to overcome some 
of the shortcomings of the previous studies and provide a reliable and a de- 
tailed description of the effects of 3d transition-metal impurities in MgCNi 3 . 

In the present study of the substitutional^ disordered MgC(Nii_ x T x ) 3 al- 
loys, with T=Fe, Co or Cu and < x < 1, the effects of chemical dis- 
order are taken into account by employing the first-principles approach of 
Korringa-Kohn-Rostoker in the atomic sphere approximation in conjunction 
with coherent-potential approximation (KKR-ASA CPA) [17]. The effects of 
incipient magnetism in MgC(Nii_ x T x ) 3 alloys are studied phenomenologi- 
cally using Ginzburg-Landau coefficients obtained from total energy versus 
magnetic moment calculations carried out with fixed-spin moment method. 

To obtain the Ginzburg-Landau coefficients one first calculates, using the 
fixed-spin moment method, the magnetic energy AE(M) = E(M) — E(0), 
where E(0) and E(M) are the total energies for the paramagnetic phase 
and ferromagnetic phase with a net magnetic moment M respectively, for 
MgC(Nii_ x T x ) 3 alloys. Then the calculated AE(M) as a function of M is 

fitted to the form ^ — a 2n M 2 ™ for n = 4. The fit coefficients a 2n can be 

n>0 2n 

used to study the propensity of magnetism in the substitutionally disordered 
MgC(Nii_ x T x ) 3 alloys as a function of x. For example, in the phenomenolog- 
ical approach the coefficient a 2 < would represent a ferromagnetic ground 
state while a 2 > would imply a paramagnetic state. 

Before describing our work in detail, we briefly describe the main results of the 
present study. We find that the disordered MgC(Nii_ x T x ) 3 alloys have a small 
magnetic moment localized at Fe and Co sites for low concentrations. The 
overestimation of the calculated magnetic moment is likely to be due to the 
limitations of the local-density approximation. From the calculated Ginzburg- 
Landau coefficients we conclude that the disordered MgC(Ni 1 _ x T x ) 3 alloys 
remain paramagnetic throughout the concentration range. At expanded vol- 
umes, we also find the possibility of a ferromagnetic state for MgC(Ni 95 Fe ,Q 5 ) 3 
and MgC(Ni s M Co , l0 ) 3 . 

The unpolarized, spin polarized, and fixed-spin moment calculations for MgC(Nii_ x T x ) 3 
alloys are carried out, self-consistently, using the KKR-ASA-CPA method. 
The calculations are scalar-relativistic with the partial waves expanded up to 
Imax = 3 inside the atomic spheres. The exchange-correlation effects are taken 
into account using the Perdew and Wang parametrization [18]. The core states 
have been recalculated after each iteration. To improve the ASA calculations, 
we have incorporated both the muffin-tin correction for the Madelung energy 
[19,20] and the multipole moment correction to the Madelung potential and 
energy [21] in our calculations. For this purpose the multipole moments of 
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the electron density have been determined up to l^ ax = 6. Further, screening 
constants a and (3 were included in our calculations following the prescription 
of Ruban and Skriver [22,23]. These values were estimated from the order(N)- 
locally self-consistent Green's function method [24] and were determined to 
be 0.83 and 1.18, respectively. The atomic sphere radii of Mg, C and Ni/Co 
were kept as 1.404, 0.747, and 0.849 of the Wigner-Seitz radius, respectively. 
The overlap volume resulting from the blow up of the muffin-tin spheres was 
less than 15%. The integration of the Green's function over the energy to 
evaluate moments of the density of electronic states was carried out along a 
semi circular contour comprising of 20 points in the complex plane. For the 
Brillouin zone integrations, the special k— point technique was employed with 
1771 k— points spread in the irreducible wedge of the cubic Brillouin zone. 

For MgCNi 3 the equilibrium state is found to be paramagnetic, in agreement 
with the previous reports [1,6,5,4]. The calculated equilibrium lattice constant 
and bulk modulus for MgCNi 3 were found to be 7.139 a. u. and 0.42 Mbar re- 
spectively, in agreement with earlier theoretical calculations [1,6]. The density 
of states at Fermi energy Ep for MgCNi 3 evaluated at the calculated equi- 
librium lattice constant is 14.56 state/ Ry — atom (5.35 state/eV — cell). This 
value has been found to be sensitive to the approximations involved in the 
band structure calculations, and have been reported to be in the range of 
4.8-6.4 state/eV - cell [1,6,5,4]. 

For MgC(Nii_ x T x ) 3 alloys for 0.0 < x < 0.2, the lattice constants decreased 
linearly for both Fe and Co substitutions at the rates of —0.076 and —0.073 
a.u./at%, respectively, while for Cu substitutions the lattice constants in- 
creased at the rate of+0.164a.u./at%. 

The spin polarized calculations for disordered MgC(Nii^ x T x ) 3 alloys showed 
a small but finite local magnetic moments at both Fe and Co sites in their cor- 
responding disordered alloys. Note that the total energies of the unpolarized 
calculations are essentially degenerate with the corresponding total energies 
of the spin polarized calculations. The variation of the local magnetic mo- 
ment in MgC(Nii- x T x ) 3 alloys as a function of x is shown in Fig.l. The local 
magnetic moments found at the impurity sites in MgC(Nii- x T x ) 3 alloys are 
inconsistent with the experiments. However, one may note that the exper- 
imental characterizations on the effects of 3c? transition-metal impurities in 
disordered MgCNi 3 alloys have remained controversial. 

The appearance of local magnetic moments in disordered MgC(Nii_ x T x ) 3 
alloys may be attributed to (i) an improper choice of the individual atomic 
sphere radius, (ii) inadequate number of angular momentum states in the basis 
set for the expansion of the wave functions, (iii) neglect of local environment 
effects such as lattice relaxations and site correlations, and (iv) an incorrect 
description of the exchange-correlation potential. 
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Figure 1. The local magnetic moments at the Fe site in MgC(Nii_ x Fe x ) 3 and Co 
site in MgC{Nii_ x Co x ) 3 as a function of x, calculated as described in the text. 

We have performed our calculations for disordered MgC(Nii_ x T x ) 3 alloys 
using different atomic sphere radii but still found local magnetic moments at 
the impurity site, namely Fe and Co. No local magnetic moments were found 
at the Cu site as well as Ni site in MgC(Ni 1 _ x T x ) 3 alloys. We have checked 
that the present calculations are well-converged with respect to (a) angular 
momentum basis set, (b) number of energy points in the contour for evaluating 
the Green's function and (c) the number of k-points in the irreducible part of 
the Brillouin zone for k-space integrations. 

The effects of lattice relaxation in disordered MgC(Nii_ x T x ) 3 alloys, if any, are 
expected to be very small. For example, the neglect of the experimentally [25] 
determined deviation of around 0.05A of Ni atoms from its ideal symmetry 
positions in MgCNi 3 in earlier [1,6,5,4] as well as the present calculations 
show no magnetism in MgCNi 3 . Moreover, first- principles calculations for 
3d impurities in a 3d host suggest that the lattice relaxation effects would be 
less important [26,27] in determining alloy energetics. Usually, these effects 
are large if there is a large atomic size mismatch between the host and the 
impurity atoms. 

The use of other exchange-correlation parametrization in our calculations of 
MgC(Nii_ x T x ) 3 alloys still leads to the formation of local magnetic moments 
at the Fe and the Co sites. Thus, the possibility of inadequate description of 
the exchange-correlation effects within the local-density approximation (LDA) 
cannot be ruled out, especially when spin-fluctuation effects are important. For 
example, FeAl, Ni 3 Ga, Sr 3 Ru20j, SrRh0 3 , and NaCoiO^ are paramagnets, 
but the LDA-based first-principles calculations show a magnetic moment of 
0.70/is [28], 0.79 hb [29], 0.64 fi B [30], 0.83/i B [31], and 0.35/i B [32,33], respec- 
tively. The Corrections to these overestimated magnetic moments are, how- 
ever, incorporated on a phenomenological level using the fixed-spin moment 
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method [34]. 



We use a similar approach to study the incipient magnetism in the disordered 
MgC(Ni 1 _ x T x ) 3 alloys. Further, Singh and Mazin [1] have anticipated spin- 
fluctuations in MgCNi 3 based on the Stoner model of itinerant magnetism. 
The overestimation of magnetic moments in the ferromagnetic state of the 
alloys within the LDA can be considered as an indicator of spin-fluctuations 
in the system. However, for this to be an effective indicator, one needs to 
rule out the competing states like anti-ferromagnetism. Empirically, it has 
been suggested [35] that anti-ferromagnetism would step in if the paramag- 
netic density of states has E F in the non-bonding region. For those systems, 
where E F is in the anti-bonding region, exchange splitting of bands leading 
to ferromagnetism is more likely. In MgCNi 3 , the Ep is in the anti-bonding 
region, and hence it is closer to ferromagnetic instability rather than to an 
anti-ferromagnetic one. The absence of nesting features in the Fermi surface 
also rules out anti-ferromagnetism in MgCNi 3 [6,36\. 

We have calculated the magnetic energy AE(M) as a function of M for various 
values of x for MgC(Nii_ x T x ) 3 alloys at the equilibrium lattice constants using 
the fixed-spin moment method. The calculated results of AE(M) as a function 
of M for MgC(Nn_ x Fe x ) 3 , MgC^Ni^Co^ and MgC(Ni 1 _ x Cu x ) 3 alloys 
are shown in Fig. 2. It can be seen that the curve AE{M) is rather flat near 
M=0 for small values of x in all the alloys and could be fit very well with the 

form —a 2n M 2n for n = 4. 

n>0 2n 

The changes in the Ginzburg-Landau coefficients function of x are 

shown in Fig. 3. From Fig. 3, we find that a 2 decreases monotonically with in- 
creasing x for MgC(Nii^ x Fe x ) 3 and M gC (Ni\^ x C o x ) 3 alloys in the range 
0.0 < x < 0.2. Since g^'s are positive for the MgC(Nii^ x T x ) 3 alloys, one finds 
that the fixed-spin moment method in conjunction with Ginzburg-Landau 
coefficients predicts the ground state of these alloys to be paramagnetic, in 
agreement with the experiments. The variation in and are also important 
as negative values would indicate a metastable ground state. In the renormal- 
ization approach to include spin-fluctuation effects from the first-principles, 
the coefficients of the lower powers of M are corrected by means of higher-order 
coefficients. However, in the present case, one may expect that the trend in the 
variation of a 2 would remain more or less similar, as the variation in and 
with respect to x compliment each other. If one takes a 2 =0 to be the magnetic 
line of instability then as an increasing function of a; in MgC(Nii^ x T x ) 3 alloys, 
one finds that both Fe and Co substitutions in low concentration enhance the 
incipient magnetic traits when compared to MgCNi 3 itself. For MgCNi 3 , the 
Ginzburg-Landau coefficients a 2 and 04 were determined to be 470.5 and 

139.2 -T- , respectively. 

Ms 
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Figure 2. The magnetic energy AE(M) as a function of magnetization M 
for MgC{Ni\^ x Fe x ) 3 (upper panel), MgC{Nii_ x Co x ) 3 (middle panel ), and 
MgC{Nii- x Cu x )z (lower panel) alloys, calculated at their respective equilibrium 
lattice constants . The curves in the figure correspond to different concentrations x 
in at% as indicated. 

Finally, we describe the effects of volume on the magnetic properties of MgC(Nii_ x Fe x )z 
and MgC(Nii_ x Co x )3 alloys. To understand the magneto- volume relationship 
in MgC(Nii- x Fe x )z and MgC(Nii_ x Co x )3 alloys, we have carried out fixed- 
spin moment calculations over a range of volume for MgC(Ni 95 Fe om ) 3 and 
M gC (Ni 90 C o Q 1Q ) 3 alloys, respectively. We have chosen these concentrations 
because the calculated local magnetic moments are maximum at these con- 
centrations. 

The calculated magnetic energy as a function of magnetization for various lat- 
tice constants is shown in Fig. 4. Fig. 4 indicates the possibility of a ferromag- 
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Figure 3. The Ginzburg-Landau coefficients (a) 02 in units of 10 3 .^-, (b) 04 in units 

of 10 3 .-^-, and (c) ae in units of 10 3 .-^-, as a function of x in MgC(Nii_ x Fe x ) 3 , 

MgC(Nii- x Co x )3, and MgC{Ni\- x T x ) 3 alloys, calculated at their respective equi- 
librium lattice constants as described in the text. 
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Figure 4. The magnetic energy AE(M) as a function of magnetization M for 
MgC{Nio$§Feo. 05)3 (upper panel) and MgC{Nio.§oC 00.10)3 (lower panel ) alloys, 
calculated at various lattice constants as indicated. 

netic state for MgC(Ni .Q 5 Fe .o5)3 and MgC(Ni .QoCo .i ) 3 alloys. We also 
find that the paramagnetic state becomes stable at smaller lattice constants 
for both Fe and Co substitutions in MgCNi 3 . However, for Cu substitutions 
(not shown), no such volume dependency on magnetism was observed. As in- 
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dicated earlier, following the rigid band picture, electron dopings in the form 
of Cu substitutions in MgCNi 3 decreases the density of states at Ep thus 
making it paramagnetic. 

In conclusion, the first-principles, local-density-functional based calculations 
for substitutional^ disordered Fe and Co impurities in the Ni sub-lattice 
of MgCNiz, in low concentrations, show that incipient magnetism resides 
in these materials. The overestimation of the calculated magnetic properties 
points to the limitations of the local-density approximation. Using a phe- 
nomenological approach based on Ginzburg-Landau coefficients and the fixed- 
spin moment method, we have shown that the disordered MgC{Nii_ x T x )^ 
(T=Fe, Co or Cu) alloys remain paramagnetic. At expanded volumes, we 
also find the possibility of a ferromagnetic state for MgC(Ni , 9 5Feom)3 and 
MgC(Ni , 90 C 00.10)3 alloys. 
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